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ABSTRACT 



The application of phase plane analysis to a missile jet 
relay control system is reviewed, explained, and illustrated. 
Only a single lead network is considered under conditions 
in which aerodynamic damping is negligible. The assump- 
tion that all forces occur as step inputs is used to simplify 
analysis techniques and does not constitute a significant limi- 
tation on the application of this method for design purposes. 

To aid in system design, an analytic method for the 
solution of phase plane plots is developed. The analytic 
technique reduces the time necessary to evaluate a system 
design by minimizing the necessity for construction of phase 
plane diagrams. 

in conclusion a summary of general principles which 
may be used to optimize a missile jet relay control system 
is presented. 



PROBLEM STATUS 

This is an interim report on one phase of the problem; 
work is continuing. 



AUTHORIZATION 
NRL Problem A02- 18 



Manuscript submitted 2 September 1958 
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INTRODUCTION 



The jet relay control system has wide application in the field of missile control and 
guidance. A jet relay control system is generally composed of a gyro reference unit, a 
control amplifier, pilot and slave relays, and solenoids which actuate reaction jets. The 
operation of this system may be described by considering control in a single axis. The 
gyro output signal, which indicates missile angular deviation from a reference (gyro 
error), is applied to the control amplifier. After modification by the amplifier the signal 
is applied to a balanced- output amplifier. The output of this amplifier is biased so that 
the magnitude of the error signal must exceed a predetermined value, plus or minus, to 
cause pilot relay operation. A heavy-duty relay is slaved to the pilot relay and when the 
error exceeds the amplifier "dead space," the action of the slave relay energizes a 
solenoid which turns on a reaction jet. The action of the jet creates a correcting moment 
which rotates the vehicle in such a direction as to reduce gyro error and cause relay 
drop-out (Fig. 1). 

This type of control system is most effective and economical where relatively 
small correcting forces are necessary— usually under flight conditions wherein 
aerodynamic disturbances are small or negligible. 



Since this method of control is 
nonlinear, frequency response anal- 
ysis techniques are not applicable, 
and a system using phase plane anal- 
ysis has been generally adopted. 
This report will describe a phase 
plane method for analysis of mis- 
sile jet relay control systems and 
will develop analytic design criteria 
based on this analysis. The tech- 
nique described in this report makes 
the simplifying assumption, valid 
for design purposes, that the forces 
acting on the vehicle occur as step 
functions. The phase plane diagram 
is a plot of the velocity of the con- 
trolled variable as a function of its 
displacement. This technique is only 
valid when the differential equation 
describing the system dynamics may 
be reduced to one of the second order. 
Missile dynamics in the absence of 
aerodynamic damping meet this 
criterion. 

An example of a phase plane 
diagram is given in Fig. 2. Con- 
sider that some disturbance of short 
duration has resulted in an angular 




, ERROR 



Figure 1 - Typical jet relay control system 
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Figure 2 - Simple phase plane 



velocity 0 O of the missile. The angular deviation of the missile from the zero reference 
will increase with time and will generate on the phase plane a trajectory, the line 0, 1 and 
2. When the angular deviation of the missile has reached the end of the dead space, rep- 
resented by the intersection of the trajectory and the pull-in line, point 1, the control 
relay will be energized. The relay actuates a solenoid and jet thrust is developed. Since 
there is a time delay between relay energization and full thrust, the angular deviation of 
the missile will proceed along the constant-velocity path during this time, from point 1 to 
point 2. It will be shown later that when the reaction jet exerts thrust the vehicle trajec- 
tory on the phase plane will be parabolic. The acceleration due to the jet action will 
reverse the velocity along the parabola from point 2 to point 3. Since relays exhibit 
hysteresis, the angular deviation at which the relay is de-energized will differ from the 
pull-in line as represented by the drop-out line. 

There will be a time delay between the de- energization of the relay and the decay of 
the thrust, during which the trajectory will continue along the parabola from point 3 to 
point 4. The system then has obtained a constant angular velocity opposite in sense to 
the initial condition and the sequences will be repeated. The "limit cycle" is established 
when the phase plane pattern is exactly repeated. 

In the above simplified example the pull-in and drop-out lines were represented by 
vertical lines which indicate that the control amplifier did not modify the error signal. 
In practice, lead networks are employed in the control amplifier in order to damp out 
transient disturbances and to increase the period of hunting for purposes of fuel 
conservation. 

This report will develop the equations of the relay pull-in and drop-out lines for a 
single lead network, the general trajectory of vehicle motion on a phase plane plot, the 
relationship between the phase plane plot and time, analytical methods for determining 
the system limit cycles, and general rules for optimization in system design. 



EQUATIONS FOR VEHICLE TRAJECTORY ON THE PHASE PLANE 



From an examination of missile dynamics in the absence of aerodynamic forces, it 
can be seen that the angular deviation 6 of the missile heading with respect to the roll, 
pitch, or yaw axis is given by the following equation: 

0 = e o + <v + - . w 

and differentiating this with respect to t gives 

8 = e 0 + a K t ....... ( 2 ) 

where 

a g = the system angular acceleration, 

8 0 = the initial angular displacement, and 

. S 0 = the initial angular velocity. 

By combining these equations to eliminate time, the angular deviation may be expressed 
as a function of angular velocity: 

■ e= . e . o+ o p (3) 
■ ' g g i 

For the simplification of the analysis, all quantities should be nondimensionalized, 
since then a parabolic template may be made for analyzing the system for various 
parameter variations. As a base for nondimensionalizing all quantities having time 
dimensions, t b = 1 second is used. This defines a nondimensional time scale r - t/t b . 

To nondimensionalize the other quantities, the control acceleration, a, is used. It 
follows that angular displacement, velocity, and acceleration will then be defined as 

n = -2- n = -£- h = ^- 
«t b a ' at b" 

Hence Eq. (3) may be written in the nondimensional form: 

n = n o + -TtZ - V) + i -75 Cn - n 0 ) 2 , ( 4 ) 
which may be simplified to: 

n = n o *\-t(& ~ d o 2 ) • (5) 



This equation represents a parabola which is symmetric about the n axis and with vertex 
at 
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There are four general conditions which may govern a vehicle trajectory in the phase 
plane: 

(a) No disturbance, no correction, a g = 0. In this case the vehicle trajectory is rep- 
resented by a constant-angular-velocity line 

n = k. 

The following three conditions result in a parabolic trajectory: 

(b) No disturbance, correction applied, a g = a. The parabola will be 

n = n o t \ (n 2 - ft D 2 ) . 

(c) Disturbance, no correction, a g = = acceleration due to the disturbance. The 
parabola will be 

n = n 0 ±^ ( ri 2 -ri o V 

(d) Disturbance, correction applied, a = (a D + a) where the sign of a will be opposite 
to the sign of n so that it will oppose the vehicle angular motion. The parabola will be 

Before proceeding with examples of these conditions, however, it is necessary to derive 
equations for the relay drop-out and pull-in lines and develop expressions for determina- 
tion of time on the phase plane. 



GENERAL EQUATIONS OF RELAY DROP-OUT AND PULL-IN LINES 



(SINGLE LEAD NETWORK) ON THE PHASE PLANE 



To determine relay drop-out and pull-ih lines on a phase plane diagram it is neces- 
sary to develop an expression for the output voltage of the control amplifier as a function 
of its input voltage. This expression then will relate the controls response to vehicle 
motion. 

For most relay control systems the control amplifier may be closely approximated 
by a single lead network. The general form for the output of a single lead network in 
terms of LaPlace transforms is: 



W-i^. (6) 
('••- S ) 

where . 

E 0 (R') = output voltage, 

E.(S) = input voltage, and 
r,T = constants of lead network. 
This expression, written as a differential equation, becomes 

dE„(t) dE,(t) 
7E 0 (t) + T — g^— = E.(t) + T * J . (7) 

This is a first order linear differential equation which has a general solution as follows: 

(8) 



E 0 (t) =■ e T 



T t 



■ j [Ei(t) +T« i (t)J e T * dtl . 
The integration may be performed by parts and the equation becomes 



The above method of expansion could be continued, but for most control problems a 
sufficient degree of accuracy may be obtained by assuming that all vehicle angular 
accelerations occur in the form of step functions, i.e., E.(t) is either constant or zero. 
Thus Eq. (9) becomes 



*i(t) +|-(r-i) Ei(t) - ^ (r-i) £1(0(1 -e*?') 



(10) 



Since the input voltage is directly proportional to vehicle angular motion, i.e., 
E A (0 = K0(t), Eq. (10) may be rewritten in terms of e, vehicle angular deviation: 
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= 0R(t) = e(t) +X (r - 1) S(t) -li(y-l) 8(t) (l-V**).. (H) 

where & R (t) is that value of 6 corresponding to the relay drop-out or pull-in point when 
there is no lead circuit in the system. Equation (11) is the exact solution for the relay 
pull-in and drop-out lines for a constant angular acceleration. In most control problems 
the transient term may be considered zero, and the approximate equation for the pull-in 
and drop-out lines may be used: 

0 R (t) = 0(t) + |-(y- 1) 9(t) (y- l) 6(t) . (12) 

Evaluation of the transient term is discussed in Appendix A. Since the transient correc- 
tion is generally insignificant it will be neglected in the remainder of this report. 

The general equation for relay pull-in and drop-out lines may be nondimensionalized 
as follows: 

„ . mfi mc fi 

^ = n+ -^-^T' (13) 

where 

T T a fL 

m ~~ (7-1). c= ~' t b = 1 sec , fi = — . 

Let fi RP be the nondimensional angular position for relay pull-in determined by the dead 
space. Then the relay drop-out value fi RD is given by 

n RD = n Rp( R ). (14) 

where K is the value of the relay hysteresis. 

It can be seen that the positions of the relay pull-in and drop-out lines depend on the 
system angular acceleration; this important point will be demonstrated in an example 
later. Moreover, since the system angular acceleration is considered a constant at any 
time, the relay pull-in and drop-out lines are represented by straight lines on the phase 
plane plot, with slopes m and intercepts determined by fi E , fi, m and c. 

There are four general conditions which affect the equations for relay pull-in and 
drop-out: 

(a) No disturbance, no correction, a R = o. The equation for the relay pull-in line in 
this case will be 

± fij^p = fi + mfi . (15) 

(b) No disturbance, correction applied, a p - a (note that a is always opposite in sign 
to fi R ). The equation for the relay drop-out line will be 

± fi RD = fi + mfi ± mc . (16) 

(c) Disturbance, no correction, a g = a n (consider > o). The equation for the 
relay pull-in line will be 

I a Dl 

±fi Rp = fi+mfi+mc — . (17) 
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(d) Disturbance, correction applied, a. g = a D + a (consider a D > o). The equation for 
the relay drop-out line will be 

± n RD = n + m n ± mc ^ + ^ , (18) 
where the sign of a will be opposite to the sign of n RP . 



RELATION BETWEEN THE PHASE PLANE PLOT AND TIME 



Time is a parameter necessary in constructing the phase plane diagram, but it does 
not appear explicitly on the plot. Therefore a relationship between time and the phase 
plane parameters must be established. 

Again, four general conditions affect the relationship between time and the phase 
plane parameters. 

(a) No disturbance, no correction, a = 0. Under these conditions, the trajectory is 
represented by a constant-angular- velocity line, fi = K. Then 

An At 



therefore 



, An 



(b) No disturbance, correction applied, a ? - a. In this case 

6 2 = 6 l + a(t 2 - tj) 
= £j + a At , 

thus 

At -4- 

a 

This expression is then nondimensionalized as follows: 

-f = n 

At = Af2 t b . 

(c) Disturbance, no correction, a ? = a D . In this case 

A5 = a^At 
An a ty, 



(19) 



(20) 



(21) 



(d) Disturbance, correction applied, a = + a. By considerations similar to the 
above, 

Antja| 
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EXAMPLE PHASE PLANE ANALYSES 



The use of the equations developed for phase plane analysis can best be demonstrated 
by examples. The following examples are not intended to illustrate design, but to present 
technique. 

EXAMPLE 1. SYSTEM WITH INITIAL DISTURBANCE ONLY 

Consider the case where a control system is subjected to an initial disturbance which 
lasts for a finite time; after this time no disturbance is present (Fig. 3). The problem is 
to determine the maximum excursion experienced by the system due to the disturbance, 
and the system efficiency in the final limit cycle. 

The initial disturbance is expressed in terms of a disturbing acceleration. The sys- 
tem is considered at rest at n = rt"= 0; and n^, n RD , a, m, c, and the system time delays 
are known quantities. 

Step 1 . Plot on the phase plane the puU-in and drop-out lines. 

(a) Since the initial pull-in may occur with a disturbance and no correction, 



(b) Until the plot has been started it will not be known whether or not the disturbance 
will still exist at drop-out. Plot the drop-out lines for no disturbance. 



(Line C) 




(23) 




o to ,:n=-±-^ft 2 , Av^n, 



I TO z: ft = ft, CONSTANT, At= K 1 

2T0 3:a=ii ( -J-(a ! -!i\) fll 

3 to 4: n = n 2 -i (n J -ft e 2 )At =ft 4 -n, 

4 TO 5. A ft, = CONSTANT 



B 



PULL-IN 



PULL-IN 
DROP-OUT 



PULL-IN 



Figure 3 - Phase plane with initial disturbance only 



9 



10 



NAVAL RESEARCH LABORATORY 



(Lines A) ±n RD = n + m n ± mc . (24) 

(c) Plot the pull-in lines needed for the limit cycle, the case where no disturbance 
exists. 

(Line B) ±n RP =n+mri. (25) 

Step 2 . Determine the length of time the initial disturbance will exist on the phase 
plane and plot the initial trajectory. 

(a) The initial trajectory will be the parabola 

fl= i§" 2 \ r; (26) 

(b) This parabola will exist for the duration of the disturbance if the pull-in line 
(line C) is not intersected during that time. The time duration A t D of the disturbance is 
known; the parabola will extend to a velocity defined by the equation 



At„ = -2- a 



D 



a- (27) 



From Fig. 3 it may be seen that the pull-in line, Line C, is not intersected during the 
time of the disturbance, point 0 to point 1. When the disturbance stops the pull-iri line 
shifts and becomes line B since there is no acceleration. In this case pull-in occurs at 
the instant the disturbance stops, since the trajectory already has intersected line B; i.e., 
pull-in is not initiated from line B but at the end of the disturbance. 

Step 3 . Because of the time delay in relay action, solenoid action, and jet thrust 
buildup, the control acceleration will not occur for some finite time after the initiation 
of pull-in, when the disturbance stops. This time delay At, is known for the system; 
during this time delay vehicle angular motion is along a constant-velocity path, point 1 to 
point 2. The length of this path is determined by the equation 

n 2 - n t = An = n x At . (28) 

Step 4 . At point 2, the vehicle angular motion is due to the control angular accelera- 
tion a, and the trajectory in the phase plane may be expressed as 

n = n 2 - \ (h 2 - n 2 2 ) . (29) 
In practice this portion of the trajectory is drawn using a parabolic template 

h = -§- n 2 . (30) 

This template is moved with vertex on the n axis until the curve intersects the trajectory 
at point 2. The trajectory is then drawn along the constant-acceleration path, point 2 to 
point 3. 

Step 5 . The maximum vehicle angular excursion is the intersection of the parabola, 
from point 2 to point 3, with the fi axis: 

^=V. (31) 



since t b = 1 second. 
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Step 6 . Point 3 is the intersection of the trajectory and the drop-out line. Again, 
during the time delay in the relay, solenoid, and thrust decay, the trajectory in the phase 
plane will continue along the parabolic path. The point at which thrust is zero, point 4, 
is determined from the equation 

At = n 4 - n 3 , (32) 

where At is the time delay, and since t b = 1 second. This procedure is then continued 
until the pattern repeats itself, at which time the limit cycle has been reached. 

With the limit cycle a, b, c, and d established, it is possible to determine the system 
efficiency. A convenient way of expressing the efficiency is by the ratio, jet-on time per 
cycle divided by the period of cycle. Hereafter this ratio will be used as the definition of 
efficiency. It should be carefully noted that the higher this ratio, called the "Efficiency 
Number," the lower the efficiency of the system: 

jet-on time per cycle /0 ,\ 

eff. no. = : — -j — j : . (o6) 

period of cycle 

With the efficiency established, the total jet-on time for any flight duration may be deter- 
mined by multiplying the efficiency by the flight time. 

Step 7 . Determine the jet -on time by taking the value of n o and f> o at a, b, c, or d. In 
the limit, the pattern should be symmetric and the absolute value of these quantities 
should be the same at any of the four points. In this case 

jet -on time = 2(2 fi 0 ) , (34) 
jet-off time = 2 ( 2 ^-)- ( 35 ) 

EXAMPLE 2. SYSTEM WITH A CONTINUOUS DISTURBANCE 

Consider the case in which a system is subjected to a continuous disturbance. If the 
acceleration due to the disturbance is large enough, the trajectory of the vehicle in the 
phase plane will oscillate about the positive pull-in and drop-out lines, as shown in Fig. 4. 
There is, however, a possibility that the vehicle trajectory will intersect the negative 
pull-in line, resulting in a phase plane unlike that of Fig. 4. Because of this possibility 
the negative pull-in and drop-out lines should be plotted. The problem, as in Example 1, 
is to determine the maximum excursion experienced by the system due to the disturbance, 
and the system efficiency in the final limit cycle. 

Again the disturbance is expressed in terms of a disturbing acceleration. The sys- 
tem is considered at rest at n = ri = 0; and n^,, n^, a, m, c, and the system time delays 
are known quantities. The phase plane for this condition is plotted as follows: 

Step 1 . Plot the drop-out and pull-in lines: 

= n + mn ± men. (36) 

(a) For the pull-in line A, Cl is dependent on the acceleration of the disturbance: 

ft lt . = X (37) 

pul 1 -in a • v ' 
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i T02: 5- n',At = Asi- ! - 



3 to 4: n, = n 




DROP-OUT 

PULL- iN x 

Figure 4 - Phase plane with continuous disturbance 



Therefore 



± ft 



RP 



ft + mft - mc . 



(38) 



Both pull-in lines should be plotted to determine the possibility of the vehicle trajectory 
intersecting the negative pull-in line. 

(b) For the drop-out line B, ft depends on the acceleration due to the disturbance plus 
the correcting moment acceleration: 



ft 



drop-out 



(39) 



where a D > 0 and the sign of a is such as to oppose vehicle motion, i.e., opposite to the 
sign of ft Kr . Therefore 



± ft RD = ft + mft + mc 



i a D + a 



(40) 



Step 2. Determine the initial trajectory due to the disturbing acceleration, 
(a) The initial trajectory will be the parabola 



2 a D 



(41) 



(b) This trajectory will continue after intersection with the pull-in line, point 1, 
until point 2. This length is determined by 



Aft = At 



(42) 



where At is the time delay of the system resulting from relay action, solenoid action, and 
thrust build-up. 
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Step 3 . At point 2 the vehicle motion will be the result of the disturbance and an 
opposing control acceleration: 

ft 2 = n i + I^ ( "2 -"i 2 >- (43) 

This trajectory will continue after intersection with the drop-out line at, point 3 until 
point 4, because of the time delay of the system: 



An = At 



a D + a I 



(44) 



Step 4 . At point 4 the vehicle motion is again the result of the disturbance. The cycle 
is now repeated, using the same parabolic equations, until a limit cycle a, b is reached. 
The maximum excursion fi m will occur at the extreme intersection of the Or axis and the 
parabola representing the disturbance plus the correction. The efficiency number is the 
same as that defined in Example 1: 



In the present case 



jet-on time per cycle IA - S 

eff. no. - !■ — — -f fL . (45) 

period of cycle 



jet-on time = 2n linit cycle |_*_| ; ( 46 j 
jet-off time-= 2fi limit cyeU \-f-\ • (47) 



ANALYTIC SOLUTIONS FOR PHASE PLANE PROBLEMS 



It is possible to determine most of the necessary design criteria for a jet- relay con- 
trol system analytically. An analytical approach to the solution of phase plane diagrams 
makes it possible to use digital computer facilities and reduces the time required to 
evaluate system design variations. 

In the following pages, equations which may be used to solve for the necessary design 
parameters of a jet relay control system will be developed and the limitations of this type 
of solution will be explained. 

SYSTEM WITH NO DISTURBANCE 

The phase plane diagram for a system with no disturbances will be of the form shown 
in Fig. 5. This is a diagram showing the system limit cycle. At point 0, jet control 
thrust becomes effective, resulting in a parabolic vehicle trajectory on the phase plane 
which intersects the relay drop-out line at point 1. From 1 to 2, the vehicle continues on 
the parabolic path for a distance corresponding to the drop-out time delay. From point 2 
to point (-3) the vehicle has a constant velocity since no acceleration exists. At point (-3) 
the vehicle trajectory intersects the relay pull-in line. Vehicle motion continues along 
the constant velocity line from point (-3) to point (-0), a distance corresponding to the 
pull-in time delay. From point (-0) where control thrust again becomes effective, the 
vehicle motion will continue in a manner similar to that previously explained and return 
to point 0. It may be shown analytically that the phase plane diagram of the limit cycle 
for a system with no disturbances will be symmetric about the and fi axes. 



n 




a 



Figure 5 - Phase plane with no disturbance 
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Alteration of the pattern of vehicle motion on the phase plane may be accomplished 
by varying the lead circuit parameters and thus changing the pull-in and drop-out line 
positions, by changing the drop-out or pull-in time delays, or by changing the control 
system dead spot. It may be possible by varying these parameters to over-compensate 
the system, in which case the system will become semicontinuous. 

The case of an over-compensated system is illustrated in Fig. 6. At point 0, jet 
control thrust becomes effective and results in a parabolic vehicle trajectory on the 
phase plane which intersects the relay drop-out line at point 1. The distance from point 1 
to point 3 on the parabola corresponds to the drop-out time delay. In this case the vehi- 
cle trajectory intersects the relay pull-in line, point 2, before completion of the relay 
drop-out time delay. Under this condition the pull-in time delay begins at point 2 and 
extends to point (-0). In other words, part of the pull-in time delay is occurring during 
the drop-out delay. At point 3 the drop-out time delay ends and the vehicle moves with a 
constant velocity from point 3 to point (-0), a distance corresponding to the remainder of 
the pull-in time delay. The pattern f rom point (-0) to point 0 is symmetric with the 
previous pattern. In the development of an analytic solution for the limit cycle of a sys- 
tem with no disturbances the over-compensated system is not considered, since this type 
of system is inefficient. 

Since the phase plane pattern for a limit cycle (Fig. 5) is symmetric with respect to 
the n and fi axes, it is possible to completely describe the properties of the system by 
knowing the coordinates of only one of the points where a slope discontinuity exists. 
Therefore, the problem of solving analytically for the properties^ a system in the limit 
cycle, reduces to the problem of solving for the coordinates ft 0> n o at point 0. It is known 
from the symmetry of the pattern that 



n 0 = n 2 = 



■ n 



n 



-2 



(48) 



4." 



(49) 



The first step is to solve for the intersection of the relay drop-out line and the parabolic 
trajectory. The equation for the relay drop-out line will be 



PULL-IN WITH ACCEL. 
PULL-IN 
DROP-OUT 




Figure 6 - Over-compensated system 
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n RD ■= fl x + mHj + mc , 
and the equation for the parabolic trajectory will be 

Solving these equations simultaneously yields 



n x = m ± 

where the minus sign applies in the present case. Now, 

+ |n 2 - Hjl = At D , 
where At D is the time delay at drop-out; since by Eq. (49) ri 2 = -fi 0 



and 



+ m =^/ii^ 



Squaring both sides of this equation and combining terms yields 

At D 2 . 

n o = —T + n o (m - At D ) - At^m + 0^ - mc 
The time delay at pull-in is 



At, 



and, at the intersection of the constant-velocity trajectory and the pull-in line, 

n 3 = n RP -mn 0 ; 

combining these two equations yields 

This expression for n o may be equated to Eq. (56) to yield a solution for n o : 

n o [2m - (Atp + Atp)] = + At^n + (n Rp - n RD ) + mc . 

For the case where the time delays are equal, i.e., At D = At p , 

At 2 



and 



n Q 2(m - At) = - ~ + At m + (fi^ - n RD ) + mc 
fi o = < n RP " mf2 0 + At p n o ) . 
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When the coordinates fi 0 , fi Q of point 0 are known, the maximum excursion in the limit 
cycle may be determined from the parabolic relationship 

n m = -% + 1 ^ • (63) 

The efficiency number of the system may be determined from n o and ri o as follows: 

jet-on time = 2(2 fl 0 ) , (64) 

jet-off time = 2 ^ j , (65) 

4 ("o + ^). (66) 



period = 4 



eff. no. = 



1 



(67) 



It should be reiterated that the above equations apply to a control system with no 
disturbances present, and that the equations are not valid for an over-compensated system. 
One indication from the analytic solution that a system has been over -compensated is the 
appearance of a negative n o . It is suggested that whenever the analytic solution is in 
doubt a quick graphical check be made. A quick means for a graphic check is to plot n o 
and ft 0 at all four points, connect the trajectories and then measure the delay times. The 
measured times should be equivalent to the known times. 



SYSTEM WITH A CONTINUOUS DISTURBANCE 



The phase plane diagram for a system with a continuous disturbance will be of the 
form shown in Figure 7. This is a diagram showing the system limit cycle. At point 0 
the control thrust becomes effective and the result is a parabolic trajectory for which the 
resultant acceleration is due to the difference between the disturbance and the control 
force. The vehicle trajectory intersects the drop-out line at point 1. From there to 
point 2 the vehicle motion continues on the same parabolic path for a distance corre- 
sponding to the drop-out time delay. From point 2 to point 3, the vehicle trajectory in 
the phase plane is again parabolic with an acceleration due to the disturbing force only. 
At point 3, the trajectory intersects the relay pull-in line. The vehicle motion continues 
on the parabolic trajectory due to the disturbance from point 3 to point 0, a distance 
which corresponds to the relay pull-in time delay. The pattern is then repeated. 

In the case of a continuous disturbance the phase plane pattern in the limit cycle is 
symmetric with respect to the fi axis. Because of this symmetry it is possible to 
describe the properties of the system completely when the coordinates at the intersec- 
tion of the two parabolas are known. 

Because of the symmetry of the pattern (Fig. 7) 

n o = - n 2 , (68) 

n o = O, . (69) 

The first step is to solve for the intersection of the trajectory with no correcting force 
applied, and the pull-in line. The equation for the pull-in line will be 
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Figure 7 - Resultant phase plane for a system with 
a continuous disturbance 



• Kl 
%> = n 3 + mn 3 - j-^j mc < 

and the equation for the parabolic trajectory will be 



Solving these equations simultaneously gives 



n 3 = - n m. 



1*1 \ 2 2|aJ "2 / ' | a 



0, 2 1 + -r-T- I ftRp + mc 



lal 



The next step is to solve for the intersection of the trajectory with correcting force 
applied and the drop-out line. The equation for the relay drop-out line will be 



n RD = rij + m n t + mc — - — 



and the equation for the parabolic trajectory in this case will be 



2 a + 



2|a + aJ 1 



Solving these equations simultaneously gives 



|a + Ojjlr 



a+ a D |m\ 2|a + a^j 



|a+ Ojj| 



a. + 



2 a+ < 



Next the time delays are considered. The drop-out time delay will be 
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n 2 - = - At D 



I a + .aJ 



By substitution from Eq. (68) this becomes 

la + aj 

The pull-in time delay will be 



(75) 



(77) 



n - At, 



= n, 



(78) 



Now considering Eqs. (68) and (69), Eqs. (72) and (78) may be combined and rewritten 



a. = 



At r 



+ mAt, 



+ ( n R P + mc |ir|) " "o< m - 



At p ) ; 



and similarly Eqs. (75) and (77): 



a. 



- — ; -mAtr, — 



+ ( n RD - mc 



2 | B | "™ n |.| V" KU "" W 

Finally, Eqs. (79) and (80) may be combined to yield a solution for n, 



+ n o (m - At D ) 



(79) 



(80) 



At,? | a + aJ Atp 2 
n o [2m - (At D + At p )l = - -f . , - 

| a| * 



5| + » I At 



a + a r 



(81) 



For the case where the time delays are equal, i.e., A D = At p , the value of n o will be 
identical to the value obtained for a system with no disturbance, 

n 0 [2(m - At)] = - + mAt + (Ogp - n RD ) + mc . 

With n o known, n may be obtained from Eq. (79): 



0. 



2 I a 



+ mAt p |-^| + (f^p + mc 1 -n 0 (m - At p ) 



(82) 



(83) 



The maximum excursion in the limit cycle may be determined by using the parabolic 
relationship 



|a| 



° 2|a + aj 



(84) 



The efficiency number as defined earlier, for a system with a continuous disturbance, 
depends only on the disturbing and control accelerations. This may be shown as follows: 



e f f . no . = 



jet-on time per cycle 



period of cycle 
period of cycle = jet-on time + jet-off time 

= 2fl_ 



(85) 
(86) 
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Since a and will always be opposite in sign the equation may be rewritten without the 
absolute value signs as follows: 



period of cycle = 2% ^ tt - Op + ^) 

(87) 



2«o 



( 2? V 



2n D ° 



eff.no. = , ^ ~y v = X (88) 



2 n / a? \ a 



SUMMARY OF EQUATIONS FOR THE ANALYTIC SOLUTION OF 



PHASE PLANE PROBLEMS 



SYSTEM WITH NO DISTURBANCE 

For a system with no disturbance the characteristic coordinates, derived originally 
as Eqs. (59) and (60), are 

At 2 

n o [2m - (At D + At p )] = - — + At D m + (f^p - Cl^) + mc , (89) 

n o = \p + *V At P - m > '• ( 9 °) 
the efficiency number of the system, derived originally as Eq. (67), is 

„ 1 
erf. no. = — ; 

1 + -r^ (91) 

A, 2 

and the maximum excursion in the limit cycle is 

n m =n 0+ lrt 0 2 . (92) 



SYSTEM WITH A CONTINUOUS DISTURBANCE 

For a system with a continuous disturbance the characteristic coordinates are 



n o [2m - (At j, + Atp)] 



At* |u + a D | Atp 2 lau 



2 I a 



At, 



+ m At 



+ mAtr 



the efficiency number of the system is 



+ fipp + mc 



n Q (m - Atp) 



eff. no. = — ; 



and the maximum excursion in the limit cycle is 



2|a. + a^l » 



(93) 



(94) 



(95) 



(96) 
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GENERAL PRINCIPLES FOR OPTIMIZATION IN SYSTEM DESIGN 



The design of a relay jet control system consists of determination of the control 
thrust and amplifier dead spot necessary to satisfy the system limitations and determina- 
tion of the lead network parameters which will result in the greatest system efficiency 
(the lowest efficiency number). Since every system will have some limiting angular 
excursion, the choice of the amplifier dead spot and control acceleration will be based on 
this limitation. 

Consider a very simple system with a disturbance, no lead network, no time delays, 
and no relay hysteresis (Fig. 8). In this case the vehicle trajectory will immediately be 
in a limit cycle and by solution of the equations for the trajectory, the system maximum 
excursion may be expressed as follows: 

where 

n D = amplifier dead spot, 
a n = disturbing acceleration, 
a = control acceleration, and 
n = maximum excursion. 

The consideration of time delays and relay hysteresis would increase the maximum 
excursion, while the use of a lead network would reduce it. Since these effects tend to 



a 




Figure 8 - Phase plane without lead circuit or 
hysteresis 
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compensate they are neglected in the 
initial design choice of control accel- 
eration and dead spot. The relation- 
ship between maximum excursion, con- 
trol acceleration, and dead spot for the 
simple system described above is 
shown 'in Fig. 9. 

For a system with an initial dis- 
turbance, which eventually oscillates 
in a limit cycle with no disturbance, 
the efficiency number as defined earlier 
becomes less (the efficiency becomes 
greater) with increasing a and also 
with larger dead spots. 




Figure 9 - Dead spot vs ratio of disturbing 
to control acceleration for different maxi- 
mum excursions 



For a system with a continuous 
disturbance the most efficient system 
has the largest possible control accel- 
eration, since in this case efficiency depends only on control acceleration. The optimum 
lead network will depend on the time delay and the relay hysteresis as well as the control 
acceleration and dead spot. 

After selection of a control acceleration and amplifier dead spot which satisfy the 
system limitations, it is necessary to determine the optimum lead network. 

For a system with only an initial disturbance, the lead network is optimized on the 
basis of maximum system efficiency in the limit cycle. Figure 10 shows the relationship 
between y and T, the lead circuit parameters, and efficiency number in the limit. From 
consideration of this figure the value of y chosen for the lead network should be the 
largest value which can be attained and still provide an acceptable system noise level. 
After the determination of y, a series of calculations employing the equations developed 
earlier, may be used to determine the optimum 
value of T . 

For a system with a continuous disturbance, 
the lead network is optimized on the basis of the 
smallest maximum excursion in the limit cycle. 
The relationship between y and T, and the max- 
imum excursion in the limit cycle are shown in 
Fig. 11. After the largest y consistent with 
system noise considerations has been chosen, 
T may be determined by a series of calculations 
employing the equations developed earlier. An 
interesting side light on this type of system is 
illustrated in Fig. 12. This figure compares 
the initial system excursion (the first excursion 
experienced by the system if it is initially at 
rest and the disturbance appears as a step input) 
with the maximum excursion in the limit. It 
may be seen that the optimum value of T on the 
basis of the maximum limit cycle excursion 
results in an initial excursion which has approx- 
imately the same value, indicating that the sys- 
tem is nearly critically damped for this value 
of T. 



ti: 

S 

< 

< 



-t + 20 
J- + 30 t 7 + 10.. 




2 3 
EFFICIENCY 



Figure 10 - Lead circuit parameter 
T vs efficiency number for different 
values of y 
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The system time delay has not 
yet been mentioned. Figure 13 
shows the relationship between the 
system efficiency, or maximum 
excursion in the limit, and the sys- 
tem time delay. It is obvious that 
time delays should be minimized as 
much as possible. 

Relay hysteresis generally has a 
relatively small effect on the system 
efficiency. 
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Figure 11 - Effect of different values of yon 
the relation between the lead circuit parame- 
ter T and the maximum excursion in the limit 
cycle, for a system with a constant disturbance 



Figure 12 - Relation between the lead circuit parame- 
ter T and the initial and maximum excursions in the 
limit cycle for different values of y, in a system with a 
constant disturbance 
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Figure 13 - Relation between efficiency number and system 

time delay 



Appendix A 

TRANSIENT CORRECTION FOR RELAY PULL-IN AND DROP-OUT LINES* 



Correction when a Constant Velocity Exists 

The actual relay pull-in boundary on each side of the zero displacement axis is 



± n RP = n + mn - m n e" c . ( A1 ) 



Since pull-in occurs along a constant fiaxis, n = 0. 

Assume that it is required to find a relay pull-in point in the lower half of the phase 
plane. Here, n is known and a negative quantity. Along a constant n path 



T _ _An n o ~ n p 
where ' 

fi 0 = the initial point on constant- velocity path, and 

n p = the actual relay pull-in point. 
The relay will pull in when the following conditions obtain: 

T 

fip - f^p = mO - mf2e ~ c 



(A2) 



= mfi \1 - e 
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(A3) 



.•. In 1 



1 t e , 

|n|c 



Adding and subtracting n^p/lhlc from the right-hand side of this equation and multiplying 
through by c/m yields 



JI in fi . °p~ "rp \ "p- "rp . p 0 -"k 

m \ mn /" m |n| \ m |n] 

Let 

np - njjp 



mi 



mini 



(A4) 

(A5) 
(A6) 



♦Adapted from unpublished work done by J. S. Pistner of The Martin Co. 
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then 



— In (1 - y) = y - x . 
m 



(A7) 



For any value of the lead circuit parameter c/m = (l/y- l), x will approach infinity 
as y approaches unity. Hence, if a curve of y-x = ln(i- y), in other words c/m = l or 
y - 2, is plotted in the xy plane (Fig. Al) there will be a point somewhere on this curve 
with the same value of y as would exist on a curve where c/m has a value other than unity. 
The value x D = n o - is known in each particular case. 



y-x=^ln (l-y) 




Figure Al - Comparison of transient correction 
curves (no acceleration) for c/m = 1 



It is required to find yj for the case when c/m is not unity. But y t will also be on 
the curve of y- x = ln(l - y) with a corresponding abscissa x 1 . Substituting the coordi- 
nates of these points in the equations of the two curves yields 



or 



= £ln(l - yi ), 
= ln(l - yj) ; 



■5)- 



m 

X "c X ° 



m 

C X c 



or 



1 . y - 1 
7 - 2 x y - 2 x o • 



This is the equation of a straight line with the following properties: 

1 



slope = - 



y- intercept 



y - 2 ' 
y - 1 x o ' 



x-intercept - (y - l)x Q . 



(A8) 
(A9) 
(A10) 

(All) 

(A12) 
• 

(A13) 
(A14) 
(A15) 
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This straight line may be drawn, since x Q and y are known in each particular case. The 
intersection of this line with the curve y-x = ln(l - y) then gives the required value of y. 

Since the slope of the straight line only depends on the value of y, for a particular 
lead circuit all lines corresponding to different values of x o will be parallel to each other. 
The procedure for finding the actual relay pull-in point when f^p, y, and x D - n o - r^p/ml Q| 
are given is therefore as follows: On Fig. A2, draw the straight line of Eq. (A12). Read 
off the value of y where the straight line intersects the curve y-x = ln(l - y). The actual 
relay pull-in point is then given by 

*=^r' < A15 > 

or 

Op = m|fi|y + Hjjp . (A16) 




0 0.4 0.6 1.2 1.6 2.0 2.4 2.8 3.2 3.6 4.0 4.4 4.8 



x 

Figure A2 - Transient correction curve (no acceleration) for c/m = 1 



Operations in the upper half of the phase plane, where fi is positive, can be analyzed 
in a similar manner. Here it can be shown that if 

x = ik^_fk, (A17) 
m|n| 

y = Sttj^- • (A18) 

m|n| 

the foregoing procedure for finding Hp may be used except that now 

Hp = - m|n|y . (A19) 
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For a lead network where y = 2, it is only required to calculate x 0 ; then y is directly 
read off the curve y - x = ln( 1 - y). 

Correction when a Constant Acceleration Exists 

To determine the actual relay drop-out or pull-in point when an acceleration exists 
consider the general equation for relay pull-in or drop-out 

Cl^ - n + mfl - mc n 

The steady- state drop-out line will be 

n R = n + mn - mc n . (A21) 

Let n A be the approximate relay drop-out point (steady- state), and Oj. the actual relay 
drop-out point; then 

Oj, - n A = mf2 A - mc n , (A22) 

n R ~n x =mft r -~mcfi(i-V" 5 '). (A23) 
Dividing both equations through by men gives 

Now Tj the nondimensional time to go from the beginning of the acceleration 8 to the 
drop-out line will be 

r-M^-f; (A26) 

therefore 



mcfl cfi 



1 - • cR J. (A27) 



In order to solve the equations, consider the phase plane axis to have been transformed 
so that n\ = 0. Then 

An = n 2 . (A28) 

With this transformation, substitute 
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Then Eqs. (A25) and (A27) become 

<C + 1) - x = y , (A31) 
(C + 1) --x = y + e-y , (A32) 



where 



a, 

C=~. (A33) 
men ' 



The general equation for the trajectory on the phase plane with an acceleration ft 
will be 

Dividing this equation through by men gives 

n 2 - nj ! c ft 2 2 



men 2 m c 2 n 2 ' 



(A3 5) 



Now substitute 



n 2 



mc: 



n 2 



Z ' (A36) 



c 



ft ' 



(A37) 



and the equation becomes 



where 



x-D = jfy 2 , (A38) 



A., 
mcft 



(A39) 



If these curves are plotted it may be seen that actual relay pull-in or drop-out occurs 
before the trajectory intersects the steady-state pull-in or drop-out curve (Fig. A3). 

Examination of the equation for the trajectory, Eq. (A38), indicates that, plotted on 
log- log paper, it is the equation of a family of straight lines of constant slope. Hence if 
Eqs. (A3 1) and (A32) are also plotted on log- log paper, the actual relay drop-out or pull- 
in point will be given by the intersection of the straight line (representing the trajectory) 
drawn through the approximate relay drop-out or pull-in point determined from the 
phase plane plot, and the line representing the actual relay drop-out or pull-in point. 

In summary, the procedure for finding the actual relay pull-in or drop-out point 
when an acceleration exists is as follows: 

1. On the phase plane plot draw the steady-state relay pull-in and drop-out lines. 

2. Find the nondimensional change in velocity, Aft, from the initial point on the 
parabola to the approximate relay pull-in or drop-out point as given by the intersection 
of the parabola and the steady- state relay pull-in or drop-out line. 
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Figure A3 - Relation between actual relay pull-in 
and drop-out and phase plane trajectory 

3. Divide |Af)| by ch. 

4. Enter this value |An|/cfi on the transient correction plot (Fig. A4) and find the 
corresponding point on the steady- state line. 

5. Through this point, draw the straight line representing the parabola. This is a 
straight line parallel to the straight line labelled as parabola (reference) on the plot. 

6. The intersection of this line with the actual drop-out line determines the actual 
drop-out point. Read off this value of |AH/cfi| . 

7. Multiply this value by cfi and adjust the value of |Af2| on the phase plane plot to 
obtain the actual relay pull-in or drop-out point. 




Figure A4 - Transient correction for relay drop-out or pull- 
in line when an acceleration exists 



Appendix B 
LIST OF SYMBOLS 



a = angular acceleration due to a control moment 
op = angular acceleration due to a disturbance 
<x K = system angular acceleration 
y = lead circuit parameter 
8 = vehicle angular displacement 
6 = vehicle angular velocity 
8 = vehicle angular acceleration 
8 o = initial angular displacement 
0 O = initial angular velocity 

8 R = value of 6 corresponding to the relay drop-out or pull-in point when there is no 
lead circuit in the system 

t = t/t b = nondimensional time scale 

n = nondimensional angular displacement 

n = nondimensional angular velocity 

n = nondimensional angular acceleration 

fij, = amplifier dead spot 

n m = maximum excursion in limit cycle 

n RP = nondimensional pull-in point determined by dead spot 

n RD = nondimensional drop-out point determined by dead spot 

c = T/y 
E 0 (s) = output voltage 
Ejfs) = input voltage 

m = T( r - l)/y 

t = time 

t b = 1 second (time) 

At = system time delay 
At D = system time delay at drop-out 
At p = system time delay at pull-in 
T = lead circuit parameter 

* * * 
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